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Abstract 

We consider three examples of weekly perturbed centers which do not have geomet- 
rical equivalence: a linear center, a degenerate center and a non-hamiltonian center. 
In each case the number and amplitude of the limit cycles emerging from the period 
annulus are calculated following the same strategy: we reduce of all of them to lo- 
cally equivalent perturbed integrable systems of the form: dH{x, y) + e{f{x, y)dy — 
g{x,y)dx) = 0, with H{x,y) = |(x^ + y"^). This reduction allows us to find the 
Melnikov function, M{h) = jH=h f^V ~ gd^i associated to each particular problem. 
We obtain the information on the bifurcation curves of the limit cycles by solving 
explicitly the equation M{h) = in each case. 
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1 Introduction 



Dissipation is present in almost every physical phenomenon. The only possibility to 
create self-sustained oscillations in a system is by re-injecting the energy dissipated 
into it. For instance, the pumping in a switched on laser must provide the necessary 
energy to compensate the lost energy ought to the interaction between the single 
frequency light beam and its surrounding (cavity, mirrors, etc.). In other context, 
the solar nuclear bomb must radiate to the space enough energy to maintain all the 
living cycles on the earth. In this case, many different frequencies are associated to 
the thousands of different species of animals and plants [|1|. Thus, every oscillation 
in nature is generated by a balance between input and output energies, in defini- 
tive between amplification and dissipation. This dynamical state (characterized by 
a preferred period, wave form and amplitude) is represented in phase space by a 
isolated closed curve called limit cycle 0. This isolated periodic motion can emerge 
from the perturbation of a conservative system where a continuum of periodic solu- 
tions exists. Only those closed orbits where the energetic balance after perturbation 
vanishes remain as periodic motions. Take, for instance, the harmonic oscillator 
X + X = 0. A pattern of circles filling the whole plane (x, x) are all the solutions 
of this system. A weak nonlinear perturbation proportional to the velocity, such as 
that in the Lienard equation x + ef{x)x + x = 0, destroys the geometrical aspect 
of the phase space. The system can balance the pumping and damping caused by 
the nonlinear term, ef{x)x, only on some special periodic orbits. In the van der 
Pol case, where /(x) = (1 — x^), it has been proved the existence, uniqueness and 
non-algebraicity of a limit cycle for every value of the control parameter e [Q. De- 
pending on the conditions imposed over the strength of the nonlinearity e and over 
the properties of f{x): degree, parity, etc., many results on the number and ampli- 
tude of the limit cycles for general Lienard systems are scattered in the literature 
(see, for instance, Ref. |^, |], |^, H, |^, |10] and references therein). This behavior 
is the mathematical evidence of how a small perturbation in a physical system can 
destroy the coexistence of infinitely many periodic motions and leaves only a finite 
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number of them verifying strict energetic balance conditions. 



Following this line of reasoning, the aim of this work is to exploit the topological 
similarities of different systems presenting an infinite number of periodic oscillations 
which give rise to the emergence of isolated closed orbits under small perturbations. 
This is an alternative perspective of the problem of identifying the limit cycles 
growing from centers respect to other methods existing in the literature [jlT], [1^, |13|, 
14| . In Section 2 we sketch the theory behind this type of systems and we develop 



the method to find those periodic solutions. In section 3 we apply this method to 
three different problems and calculate the number and amplitude of the limit cycles 
emerging from their perturbations. The comparison with previous investigations on 
these systems is also performed. Finally, we present our conclusions. 



2 Geometrical and Topological Equivalence 



The plane is the natural space for the representation of the integral curves of a 
second order one-dimensional or first order two-dimensional dynamical system. A 
periodic motion of this system draws a closed orbit on the plane. If the relation 
H{x, y) = h, with h a fixed real number, is verified by every point (x, y) on this type 
of orbit, a continuum set of periodic orbits is obtained when h runs on an interval 
of real values, /i G [hi, h2]. This orbit structure receives the name of period annulus. 
Suppose, without lose of generality, that the boundary orbit given hy H = hi is 
a degenerate orbit, say the origin (0,0). This equilibrium point surrounded in its 
immediate neighborhood by closed paths is called a center Obviously, every 
orbit of the period annulus verifies the differential relation 

dH = 0. (1) 

Depending on the time parameterization of the curves H{x, y) = h, different 
dynamical systems having the same pattern of integral paths are obtained. At a 
first sight, the difficulty to find the integrating factor in each case prevent us from 
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establishing the phase plane geometrical equivalence between them: we say that two 
systems are phase plane geometrically equivalent when their integral curves on the 
plane are the same. Think, for instance, on the set of circles Hc{x, y) — \{x^ -\-y^) — 
h. They obey the differential relation dHc = 2xdx + 2ydy = 0. Different time 
parameterizations of this expression produce systems with different time behaviors 
(e.g.: {1) X + X — ; {2) x — ky, y — —kx, k — etc., and (3) x — y^, y — —xy), 
although with the circles as their integral curves. If the interest does not reside on 
the time evolution of the system, the last concept of geometrical equivalence will 
take importance in the study of the system. 

One step further. It is also possible that two different planar systems (having each 
one a period annulus, given by H{x, y) — h and H{x, y) — h, with h e [hi, /i2] and 
h & [hi, h2\, respectively) which are not phase plane geometrical equivalent, are yet 
topologically equivalent. By this we mean that there exists a bijectivc and continuous 
transformation of coordinates (homomorphism) , F : R^, transforming one 

period annulus into the other one. That is, if we represent the integral curves of the 
first system as Ch — {{x,y) \ H{x,y) = h} and of the second one as = {{x,y) \ 
H[x,y) = h} then T{Cfi) = C/^ in a continuous and monotone way when h and h 
run in their intervals of existence. In particular, the boundary conditions over F are: 
^{ChJ — Cfi^ and ^{0^^ — Ch^. Take for example the system x — ay,y — —bx with 
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a,b > 0. It presents a period annulus formed by ellipses, H{x,y) = ^ + ^ = h. 
Clearly, this system is topologically equivalent to another one having a circular 
period annulus: Hc{x,y) — ^{x'^ + y^) — h. It is straightforward to verify that the 
function V{x,y) = i^,-^) establishes the homomorphism between both systems. 
Obviously, if two systems are geometrical equivalent they are topologically equivalent 
because in this case H = H and then we can choose F as the identity mapping. 

The relevance of this property is the following. Imagine that we are able to find 
the limit cycles that emerge from the perturbation of the first system, 

dH + perturbation = 0. (2) 
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Then, the apphcation of T to the second perturbed system, (iif + perturbation' = 0, 
bring this last expression to the form (H), even if we do not know H. Now the 
hmit cycles of the second system can be found by solving equation (0). Undoing the 
change of variables with T~^, we will obtain the periodic motions of the perturbed 
system in the original coordinates. 



3 Limit Cycles emerging from some Centers 

3.1 Perturbation of a circular period annulus 

As an application of the method above explained, we study here the limit cycles 
emerging from the perturbation of what we could consider a paradigmatic case in 
the center's typology: a continuum set of circles, Hc{x,y) = + y"^) = h. The 
solution of many different dynamical systems is this kind of period annulus. As a 
simple representation of them we take the linear center: 

^ " ^' (3) 
y = -X, 

whose integral curves clearly verify dHc = xx + yy = 0. This path diagram is de- 
stroyed when we slightly perturb this conservative system with two general nonlinear 
terms, f{x,y) and g{x,y), controlled by the parameter | e |<^ 1, 

X = y + ef{x,y), 
y = -x + eg{x,y). 

It is a time parameterized version of the differential relation: dH^+e (/ dy—g dx) = 0, 
in which the time variable has been removed. Following the procedure proposed in 
the previous section, we need to study this differential relation in order to find the 
limit cycles of system (^. If its integral curves are represented by y{x) and we define 
y'{x) = dy/dx, we obtain: 

{yy' + x) + e[f{x, y)y' - g{x, y)] = 0. (5) 

If we suppose the origin (0,0) is the only fixed point of equation (|), a limit cycle 
Ci = {x,y±{x)) of it, with a positive branch y+{x) > and a negative branch 



y-{x) < 0, cut the x-axis in two points (— a_, 0) and (0+, 0) with a_, a+ > 0. Every 
hmit cycle Ci solution of equation encloses the origin and the oscillation x runs 
in the interval — a_ < ,t < a_|_. 



The amplitudes of oscillation a_, 0+ identify the limit cycle. The result for e 7^ 
is a nested set of closed curves that defines the qualitative distribution of the integral 
curves in the plane {x, y). The stability of the limit cycles is alternated. For a given 
stable limit cycle, the two neighboring limit cycles, the closest one in its interior 
and the closest one in its exterior, are unstable, and conversely (see Figure 1). If we 
determine the stability of the origin, then the stability of every limit cycle remains 
fixed by the alternation property. The stability of the origin is determined by the 
sign of the real part of the eigenvalues of the Jacobian at this point. Hence, its 
calculation shows that if e(/x.(0,0) + gy{0,0)) < then the origin is stable and if 
e(/x(0, 0) + gy{0, 0)) > then it is unstable. 

Another property of a limit cycle can be derived from the fact that the mechanical 
energy E = He = + y^)/2 is conserved in a whole oscillation: 

/ dHc= —dx = 0. 
Jct Jci dx 

Thus, if equation (H) is integrated along a limit cycle, between the amplitudes of 
oscillation, we obtain: 

la^ [9{x, y+{x)) - f{x, y+{x))y'_^{x)]dx + 
Ia+[9{x,y-{x))-f{x,y_{x))y'_{x)]dx = 0. 

Every couple of solutions {y^{x),y^{x)} of equation (|^), vanishing in the extremes, 
?/+(«_) = y^{aJ) = ?/+(a+) = ?/_(a+) = 0, and verifying equation (^, constitute 
the finite set of limit cycles of equation (^). As a consequence of the continuity 
in the parameter e, these conditions are also valid for e = 0. For this parameter 
value, the limit cycles are circles: y+{x) = \J a?- — x'^ and y-{x) = —y/a"^ — x^. Here 
a_|_ = a_ = a > represents the amplitudes of the limit cycles. At this order of 
approximation, the condition (|^) is verified, and condition (|^) reads: 

ra 
J —a 



r-a 
Ja 



dx+ 

(7) 
dx = 0. 



that is, 

f{x, Va^ - x2) 



where 



^(x, V a'^ — x'^) + X- 



— x^ 



dx = 0, (8) 



/(a;, y) = + fix, -y) - f{-x, y) - f{-x, -y)], .g. 

9{x, y) = Udi^^ y) - 9{x, -y) + g{-x, y) - g{-x, -y)]. 

If (3{a) does not vanish identically, each solution a > of the equation (3{a) = is 
the amplitude of a limit cycle of the system (Q) in the weak nonlinear regime. And 
conversely, at order zero in e, the amplitudes of all the limit cycles emerging from 
the period annulus are solutions of equation (|^). These results are exact for e = 0. 
Therefore, equations or (|]) determine the amplitudes of the period motions 
surviving to a slight perturbation of a period annulus formed by circles. A typical 
portrait of limit cycles in this regime is given in Figure 1. 

We remark that /5(a) is the first Melnikov function of the system (^. Here it 
has been obtained by following a different approach to the usual one (based on the 
calculation of the displacement function of the first return mapping on the period 
annulus). Observe also that f{x, y) is an odd function of x and an even function of 
y, whereas g{x,y) is an even function of x and an odd function of y. Therefore, if 
f{x,y) and g{x,y) are polynomials in x and y, only the odd terms in x and even 
in y of f{x,y) and the even terms in x and odd in y of g{x,y) survive in (|^) and 
contribute to /5(a) in (§). 

As an example, we integrate equation (||) when f{x, y) is a polynomial of degree 
2ni or 2ni + 1 in x and degree 2mi or 2mi + 1 in ?/ and g{x, y) is a polynomial of 
degree 2n2 or 2n2 + 1 in x and degree 2m2 or 2m2 + 1 in Then, 

^(x,y) = E"ioi:rio&.>^v+'- 



where a^ ^ and hj^k are real coefficients. The result is 



^ a^^^ r(j + l/2)r(fc + l/2) 



1\ /, 1 



a2(^-+^-). 



where n =Max{ni,n2} and m =Max{mi,m2}. Here, aj^k = for j > rii or k > mi 
and bj^k = for j > n2 or k > m2- The root a = of /3(a) corresponds to the fixed 



point (0, 0) and the factor a? can be eliminated. Thus, the possible amplitudes a are 
the zeros of (3{a)/a?, a polynomial in a? of degree Max{ni + mi, + ^2}. There 
are no more than Max{ni + mi, 77-2 + m2} different solutions a > and therefore, 
the maximum number of limit cycles in this case is Maxjrii +mi, n2 + m2} as it was 



Max{deg f, deg 3} — ! 
2 



where 



proved by Iliev |T^. This result can also be written as 
means the integer part and deg f (resp. deg g) denotes the degree of /(x, y) (resp. of 
g{x,y)). If we restrict ourselves to the case of Lienard equations where f{x,y) = 
and g{x, y) is linear in y, we recover the known result that the number of limit cycles 
is less or equal than n2 in the weak nonlinear regime. The extension of this relation 
for the whole range of e is the still not proved Lins-Melo-Pugh conjecture on the 
number of limit cycles for Lienard systems 0. 



3.2 Reduction of some examples to a perturbed circular pe- 
riod annulus 



Example 1: (a particular circular period annulus) 

X = y - e{aix + a2x'^ + asx^), 
y = —X — eb-^x^, 



Yuquan and Zhujun have investigated the number of limit cycles of a similar 
and more general system. The conditions of their theorems can be translated for 
this particular case of weak perturbations (| e |^ 1) as following: (i) If eb^ > 0, 
eai < and ea^ > then system ( |To|) has a unique stable limit cycle {Theorem 1 
in Ref. fl^); (ii) If ebs > 0, eai > 0, and ea^ < then system ([lOD has a unique 
unstable limit cycle ( Theorem 2 in Ref. [1^ ) . 

These results are confirmed by the method introduced in the previous section. 
Moreover, as it is explained in that section, the amplitudes a of the limit cycles 
in the weakly nonlinear regime are the nontrivial positive solutions of the equation 
/5(a) = 0. In this case, f{x,y) = —aix — a2x'^ — a^x^ and g{x,y) = —b^x^. We 
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obtain: 



13(a) = 2 ^ + 4^3« 



which shows that the system has at most one hmit cycle of amphtude 



4ai 

Therefore, the system has not hmit cycles if aia^ > and just one limit cycle if 
aitts < 0. This periodic motion is stable if eai < and unstable if eai > 0. Numerical 
simulations are in strong agreement with this analytical result. Observe that there 
are certain conditions of Yuquan and Zhujun theorems that can be removed in this 
regime in order to have just one limit cycle. 



Example 2: (a system topologically equivalent to a circular period 
annulus) 



X = —y'^^ ^ + eP{x, y), 
y = x^'^'-i + eQ{x,y), 



where k and / are positive integers and P{x, y) and Q{x,y) are polynomials in x and 
y. The number of limit cycles of this system for | e |^ 1 has been studied in Ref. 



18, 19|. The constant of motion for the unperturbed problem is H{x,y) 



1 ^2k 



2k 



J. „,2/ 
21 



y . When fc, / > 1 it corresponds to a degenerate center topologically equivalent 



to the circular period annulus obtained when k = I = 1 (then Hc = \ (x^ + y"^))- If 
we remove the time variable, we obtain 

dy 



y''-'^ + x''-' + e 
ax 



Q{x,y) - P{x,y) 



dx 



0. 



(12) 



In order to apply the method introduced in the previous section, we perform the 
change of variables T : (x, y) {X, Y): 

X = sign{x)^-^, 
Y = signiy^-^. 



Then, equation (|T2|) reads 



X |(Vfc)-i 

A;l-{l/2fe) 



I Y KVO-1 dY 



(13) 
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with 

x = stgn{X){Vk\X\Y/^ 

y = sign{Y){^l\Y\YlK ^ ' 

Equation ([13D has the form (|^) with (x, y) replaced by (X, Y) and 

/(X,r) = -M;(^p(,,^^(x)(v^ I X \r/^stgn{Y)iVl \ Y 
giX,Y) = -S^Q {stgniX)iVk \ X szgn{Y){Vl \ Y . 

Therefore, in these coordinates, if A is the amphtude of a hmit cycle {X,Y{X)) in 
the weakly nonlinear regime, it verifies the equation l3{A) = 0, where 

p{A) ^ J^^ 1^ (X, 1^4 (X)) + (X, r^(x))| dx, 

f{X,Y) and g{X,Y) are defined in (|) with / (X, F) and ^ (X, F) given in (|15|) 



and Y^(X) = \/ A"^ — X^. If we undo the change of variable (X, F ) {x,y) in the 
above integral and write a"^^ = kA"^, then we can write l3{A) = (3{a) in terms of the 
amplitudes a of the original variables (x, y): 

(3{a) = f |q - P(x,ya(a;)) rfx, 

J-a [ ax j 

where ya{x) = ^ {cP'^ — x^^^ and P{x,y) and Q{x,y) are defined in for 
f{x,y) = P{x,y) and g{x,y) = Q{x,y). Hence, the amplitudes, a, of the limit 
cycles in the weakly nonlinear regime are the nontrivial solutions of the equation 
P{a) = 0. 



li k = I, then (3{a) reads 

( ~ x'^^~^y{x)P {ax,ay{x)) 

1 — X 



(3{a) = a I \Q (ax, ay{x)) H ^ > dx, 



where y{x) = (1 — x \ . If P{x,y) and Q{x,y) are polynomials of degree at 
most n, a~'^(3{a) is a polynomial of degree [^^J in a? and we conclude that the 
maximum number of limit cycles in this regime is [^^J , as it has been shown by 



Coll et al. 19 . 



Take, for instance, the particular case: 



X = —y^ + ebxy"^ 
ij = x^ + ecx'^y^. 



(16) 
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with k = I = 2 and n = 5. The equation f3{a) = has at most one non trivial 
solution given by 

~ r(3/4) V c ' 

In this case, the system has at most one limit cycle of amplitude a emerging from 
the slightly perturbed period annulus. This analytical result is in agreement with 
the direct verification of the dynamics by integrating the equations (|16]). 

Example 3: (a system topologically equivalent to a circular period 
annulus) 



X = —y + yx^ + eF{x, y) 
y = X + xy"^ + eG{x, y), 



(17) 



where F{x,y) = a^x^ + aix'^y + a2xy'^ + a^y'^ axidG{x,y) = h^x^ + hix'^y + h2xy'^ + hzy^ . 
It has been shown that this system has, at most, two limit cycles [theorem 2.6 in 
Ref. 0). 

The constant of motion of equation (|T^ is H{x, y) = = h. A period annulus 



topologically equivalent to He = {x^ + Z/^) = h is defined in the region — 1 < x < 1 
when h runs from 1 to oo. The origin (0, 0) is found for h = 1 and the boundaries of 
this region are the lines x = ±1 which are approached when h oo. The interior 
is formed by closed curves topologically equivalent to circles. The integral paths 
exterior to the period annulus are open curves obtained when h increases from — oo 
to zero. 

Following the same steps as in the preceding example, we proceed to remove the 

dx . 
dy ■ 



time variable of system ([T7|). Writing x' ~ '^^ 



(1 + y'^)xx + (1 - x^)y + e[G{x, y)x - F{x, y)] = 0. (18) 

The homomorphism with the circular period annulus is established by means of the 
change of variables T : (x, y) (X, Y): 



X = Sign{x)J — log(l — x 



Y = Sign{y)J\og{l + y"^). 
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Then, the equation ( ]TB| ) reads 

G (^Sign{X) y/l-e-^^ ,Sign{Y) e^^ -ij 



+ F + 6 



X 



dX 



0. 



Now this equation has the form (|]) with the variables {x, y) replaced by the variables 
(F, X) and 

G(^Sign(X)\/ l-e~^^ ,Sign(Y)\/ e'^^ -ij 



f{Y,X)=\X\ 
giY,X)=\Y\- 



l-e- 



l-e-x^ ,Sign(Y)y/ e^'^ -l] 



(19) 



Therefore, we can apply the method of the previous section to this equation. Then, 
the amplitude A of a limit cycle {Y,X{Y)) of this system in the weakly nonlinear 



regime, X{Y) ~ VW^^Y^, satisfies P{A) = 0, with 

= aox^ + a2xy'^ and G(x, y) = bix^y + b^y^. If we undo the change of variable 
(y, X) — s> (x, in the above integral and write = log(l + a^), then we can write 
P{A) = P{a) in terms of the amplitudes a of the original variables {x,y): 




After straightforward computations we obtain that the real roots of f3{a) are the 
real solutions of the equation /3(a) = 0, where 



P{a) = [h + 02(1 + a2)] [2 + a^- 2^1 + 02] + 
63 [2(1 + a2)3/2 _ - 2] + ao[2VTT^ - 2 - + 



We replace the variable a in this equation by the variable B = vT+c?. We look 
for solutions a > 0, that is, B > 1. The two only possible solutions are 



, _ ~(^o + ^3) ± v(ao + ^3)^ - (flo + a2)ibi + h) 

B = . 

ao + 02 

Writing 

U = B2±fe3 ^ ^ bi±bi (20) 
(10+012 ' ao+a2 ' ^ ^ 
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we have that = —u ± yuF—v. Then, the number of hmit cycles in the weakly 
nonlinear regime is 0, 1 or 2 depending on oq, 02, &i, 63 (see figure 2): 



(i) If u > — 1 and f + 2u + l>0ort;> n^, neither neither B~ are greater 
than 1. Therefore, the system has not limit cycles. 

(ii) If f + 2m + 1 < 0, then B^ > 1, B^ < 1. Therefore, the system has a unique 



limit cycle with amplitude cl = \J {B^Y " ^• 

(iii) If f + 2m + 1 > and v < u'^, then B~^, B~ > 1. Therefore, the system has 



two hmit cycles with amplitudes a"^ = y {B~^y — 1 and a = \J {B Y — 1. 

Computer integration of equation (^) is in agreement with these analytical 



results. Theorem 2.6 of Ref. pO| states that the maximum number of limit cycles 
of this system is two. By applying the method proposed in Section 2 we obtain a 
more detailed information: the exact number and amplitudes of the limit cycles of 
this system as a function of the parameters {0^,6^ | k = 1,2,3}. 



4 Conclusions 



Nowadays it is well known that the inclusion of a slight dissipation in a conservative 
system destroys almost completely the geometrical aspect of the phase space of the 
unperturbed motion. But the analytical methods available for the study of such 
perturbed systems suffer a lack of predictive power when the nonlinearities start to 
dominate the dynamics. 

In this work, we have carried out the calculation of the precise number and 
amplitude of the periodic orbits surviving to a weak perturbation of different topo- 
logically equivalent period annulus. First, the equation j3{a) =0 containing that 
information about the limit cycles emerging from a general slight perturbation of 
a continuum set of circles is stated. Essentially this equation corresponds to the 
first Melnikov function associated to the system, although it has been obtained by 
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an alternative line of reasoning. Second, the dynamical equations of other topo- 
logically equivalent perturbed centers have been reduced to equivalent equations of 
the former case. Then, the calculation of the number and amplitude of their limit 
cycles has also been carried out in the same way. The numerical integration of these 
different systems supports the analytical calculations and confirms the theoretical 
predictions. 

We must stress that, on the one hand, our method has given the correct maximum 
number of limit cycles in all the examples analyzed. On the other hand, it has 
allowed us to predict the exact number and amplitude of the limit cycles emerging 
from different topologically equivalent perturbed period annulus. It seems that 
in order to find the number and amplitude of limit cycles, it is not important to 
maintain the condition of a planar system to be a polynomial system. Perhaps 
it would be more interesting to try to identify the class of all its geometrical and 
topologically equivalent systems and to perform the calculations in the simplest 
system of this class. 
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Figure Captions 



Figure 1: A typical phase portrait of equation (^. Tlie limit cycles of ampli- 
tudes (—01,02), (—61,62), . . ., enclose the origin. Stable and unstable limit cycles 
alternate. The order of this alternation depends on the stability of the origin. 

Figure 2: The complete bifurcation diagram of system (|1^. The system has 
no periodic solutions in the white region, one limit cycle in the grey region and two 
limit cycles in the black region. See example 3 for details of the calculation of {u, v) 



from the original parameters of equation (17). 
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